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Summary. — Adopting the convection-difTusion model for energetic electron and 
proton propagation, and accounting for all the relevant hadronic and leptonic pro- 
cesses, the steady-state energy distributions of these particles in the starburst galax- 
ies M82 and NGC253 can be determined with a detailed numerical treatment. The 
electron distribution is directly normalized by the measured synchrotron radio emis- 
sion from the central starburst region; a commonly expected theoretical relation is 
then used to normalize the proton spectrum in this region, and a radial profile is 
assumed for the magnetic field. The resulting radiative yields of electrons and pro- 
tons are calculated: the predicted >100MeV and >100GeV fluxes are in agreement 
with the corresponding quantities measured with the orbiting Fermi telescope and 
the ground-based VERITAS and HESS Cherenkov telescopes. The cosmic-ray en- 
ergy densities in central regions of starburst galaxies, as inferred from the radio and 
7-ray measurements of (respectively) non-thermal synchrotron and 7r°-decay emis- 
sion, are Up = 0(100) eVcm"'', i.e. at least an order of magnitude larger than near 
the Galactic center and in other non-very-actively star-forming galaxies. These very 
different energy density levels reflect a similar disparity in the respective supernova 
rates in the two environments. A L-y oc SFR^'"^ relationship is then predicted, in 
agreement with preliminary observational evidence. 

PACS 95.85. Pw: 7-ray sources - . 

PACS 98 . 54 . Ep : starburst galaxies and infrared excess galaxies-. 

PACS 98 . 56 . Ne : spiral galaxies (M31 and M33) . 

PACS 98 . 56 . Si : Magellanic clouds and other irregular galaxies-. 

PACS 98.70.Sa: Cosmic rays - . 



1. Introduction 



In the nuclear regions of starburst (SB) galaxies, active star formation (SF) powers 
emission of radiation directly by supernova (SN) explosions and indirectly by SN-shock 
heating of interstellar gas and dust, as well as from radiative processes involving SNR- 
accelerated cosmic-ray electrons (CRe) and protons (CRp). 

Some basic considerations suggest that the timescales required for CRp to be accel- 
erated (by SN shocks) and lose energy (via pion decay into photons and e'^e~ pairs, or 
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via advection) are shorter than timescales of SB activity in galaxies, that are themselves 
comparable to galactic dynamical timescales. A consequence is that in a SB region a 
balance can roughly be achieved between energy gains and losses for galactic CRs during 
a typical burst of SF [1]. Under basic hydrostatic and virial equilibrium conditions in 
a galaxy, a minimum- energy configuration of the field and the CRs may be attained. 
This implies that energy densities of particles and magnetic fields can be in approximate 
equipartition [2]. 

The equipartition assumption enables deduction of the CRp energy density. Up, from 
the measured synchrotron radio emission (which can be observed relatively easily) and 
a theoretically motivated injection p/e ratio. Alternatively, Up can be estimated also 
from SN rates and the fraction of SN energy that is channeled into particle accelera- 
tion. Knowledge of Up enables prediction of 7-ray emission [cither at high energies (HE: 
>100MeV) or at very high energies (VHE: >100GeV)], which is mostly due to CRp 
interactions with ambient gas protons, via tt*^ decay. 

2. HE emission from star- forming galsixies 

In this section we will review some basic features of SB modeling, notably applied 
to the local galaxies M82 and NGC253, and the status of observations of star-forming 
galaxies in the HE/VHE 7-ray domain. 

2'1. Modeling. - In both nearby SB galaxies, M 82 and NGC 253, the central SB region 
(which will also referred to as the source region) with a radius of '^300 pc and height of 
300 pc is identified as the main site of particle acceleration. Here, the injection particle 
spectrum is assumed to have a non-relativistic strong-shock index q = 2. A theoretical 
Np/Nc ratio, predicted from charge neutrality of the injected CRs, is likely to hold in 
this source region - as is also the assumption of equipartition. 

Due to the implicit dependences in the expression for the synchrotron flux, CR/field 
equipartition is implemented iteratively to solve for (primaries plus secondaries), A^p, 
and B. For both M82 and NGC 253, central values Bq ~ 200 have been obtained in 
detailed models (M82: [3, 4]; NGC 253: [5, 6, 7]). 

Adopting the convection-diffusion model for energetic electron and proton propaga- 
tion, and accounting for all the relevant hadronic and leptonic processes, the steady-state 
energy distributions of these particles in the galaxies M 82 and NGC 253, in both the SB 
nucleus and the disk, can be determined with a detailed numerical treatment. In par- 
ticular, to numerically follow particle energy losses and propagation outside the source 
region, one needs to know of the HI and HII densities and their profiles, as well as the 
spatial variation of the mean strength of the magnetic field. (E.g., assuming magnetic 
flux freezing in the ionized gas, then B oc n^'^jj [8].) 

A measured radio index of ~0.7 in the central disk implies q ~ 2.4 there. This implies 
a substantial steepening of the CRe spectrum from the injection value, q — 2. The 
steady-state electron and proton spectra in the SB region of NGC 253 are shown in fig. 1. 
At low energies both spectra are flat, whereas at i? >> 1 GeV the stronger electron 
losses result in steeper electron (than proton) spectra. 

Electron emissions by bremsstrahlung and Compton scattering are shown in fig. 2; 
also shown is 7-ray emission from tt*^ decay (following pp collisions). As expected, the 
losses due to bremsstrahlung dominate the lower energy regime, whereas losses due to 
TT*^ decay dominate at higher energies. (While synchrotron emission extends to the X-ray 
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Fig. 1. - Properties of the emitting particles in the central SB region of NGC 253 [7]: steady-state 
spectra of primary (solid line) and secondary (dot-dashed line) electrons and protons (dashed 
line). 
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Fig. 2. - Properties of the emitted radiation in the central SB region of NGC 253 [7]: Radia- 
tive yields from electron Compton scattering off the FIR radiation field (dotted line), electron 
bremsstrahlung off ambient protons (dashed line), 7r° decay following pp collisions (dashed- 
dotted line), and their sum (solid line). 
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region, it is negligible at much higher energies.) Our main interest here is the VHE 
emission, which - as anticipated - is mainly from the latter process. 

It is interesting to note that the numerically predicted VHE 7-ray flux is a factor ^6 
lower than that obtained in an approximate treatment where the impact of radial energy 
losses and propagation mode of the CRp is ignored [3] . This simplified approach results 
in an unrealistically high contribution of the main disk (exterior to the SB region) to the 
total TeV emission: in fact even the slight steepening of the CRp spectrum, occurring in 
the disk because of energy losses, significantly lowers the TeV emissivity [9] and hence 
the predicted source flux. 

We note that the related neutrino flux {tt^ eventually produce + + i^e + t'/x 
+ Dfj, ) at energies higher than 100 GeV is about a third of the corresponding photon 
flux [3]. 

2'2. j-ray detections. ~ The two local SB galaxies M82 and NGC253 are the only 
non-AGN extragalactic sources that, up to now, have been detected in both the GeV [10] 
and TeV [11, 12] domains. The measured fiuxcs and spectra of both galaxies in the two 
bands agree with predictions of recent numerical models (M 82: [3, 4]; NGC 253: [5, 6, 7]). 
The highcst-SFR galaxy in the nearby universe, Arp 220, was undetected by MAGIC [13]. 

HE 7-ray detections were obtained for a number of low SFR galaxies: the Large Mag- 
ellanic Cloud (LMC: [14]), the SmaU Magehanic Cloud (SMC: [15]), the Andromeda 
galaxy (M31: [16]), and the composite Sy2/SB galaxies NGC 1068 and NGC 4945 [17]. 
The scenario of mostly hadronic HE 7-ray emission is generally confirmed for these galax- 
ies, except for NGC 1068 where emission from the active nucleus may be dominant. Only 
flux upper limits exist for the Triangulum galaxy (M33: [16]). 

3. CRs in star-forming galaxies 

The CRp energy density in galaxies can be either 1) measured directly if the GeV- 
TeV spectral flux is known, or 2) evaluated indirectly if source size, distance, and radio 
spectral index and flux are known, and particles/field equipartition and a p/e ratio arc 
assumed, or 3) estimated if the SN rate, the CRp residence timescale, the energy per SN 
going into CRs, and the size of the SF region are known. 

1) CRs and GeV-TeV emission. The detection of M82 and NGC 253 confirmed values 
Up — 0(100) eVcm"'^, resulting from accurate numerical treatments based on the solu- 
tion of the diffusion-loss equation for the accelerated particles (see refs. above). In the 
^200 pc region of the Galactic center Up ~ 5eVcm~^ based on HESS observations [18]). 
For the comparable environment of the Andromeda galaxy, the Fermi/LAT detection 
implies an average Up at a level of « 0.35 the average Galactic value [16]. For the even 
more quiet environments of the LMC and the SMC), actual Fermi/LAT GeV detections 
imply, respectively. Up - (0.2 - 0.3) and - O.leVcm-^ [14^ 15]. 

2) CRs and radio emission. Based on method 2) above, for the central SB regions of 
NGC 253, M82, and Arp 220, respectively, [1] evaluate Up « 75, 97, and 520eVcm-3. 
These radio-based estimates match, to within a factor of w2, those derived from GeV/TeV- 
based measurements. 

3) CRs and supernova rates. The CRp residence timescale is given by t"^ = t^^ +''"outJ 
where the pp interaction timescale Tpp is a function of the ambient gas density n, and 
the advection timescale Tout is a function of the speed of the outfiowing gas (fout) and 
of the size (radius) of the SB region (rg). Typically, Tpp ^ 10^ yr and ~ 10^ yr for the 
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central SB regions of Arp220 and, respectively, the M82 and NGC253 [1]. If however a 
fast SB-driven wind advects the energetic particles out of the disk plane, then possibly 
Tout << Tpp- For M82, Wout ~ 2500 kms"^ [19]. Assuming a homogeneous distribution of 
SNe within the SB nucleus of radius rge, the outflow timescale is then Tout ~ 3 x 10^ yr. 
So in some SB galaxies Tj-os Tout- During Tres a number i^sN^Ves of SN explode and 
deposit the kinetic energy of their ejecta, E^j = 10^^ erg, into the interstellar medium. 
The Galactic CR energy budget and SN statistics suggest that rj ^ 0.05 of this energy 
may go into accelerating particles. The CRp energy density in the central SB region is 
then Up = ^i'sN''"ros'7£'cjTgB . A substantial agreement with the equipartition estimates is 
reached for NGC 253, M 82, Arp 220, Milky Way, LMC: Up ~ 75, 95, 505, 5.7, 0.2eV cm'^, 
respectively [1]. Summarizing, the CR energy densities estimated in five galactic nuclei 
of similar size (three SB galaxies, the central Galactic region, and the LMC), appear to 
be largely correlated with key features of the ongoing SF: SN rate and CRp residence 
time (the latter being, in turn, a function of the local gas density and of the galactic 
superwind speed). In fact, in these environments, for same rjEcj (by assumption) and 
similar (from observations), the CRp energy density seems to be well described just 
as a function of the number of SN explosions during the CRp residence timescale. 

Up OC Z^SNTres. (1) 

4. HE/VHE7-ray emission and SFR 

The Schmidt-Kennicutt (SK) law of SF, Esfr oc (where gas comprises both HI 
and H2), states that projected SFR varies as a power law the projected gas density. This 
is true piecewise: both locally and disk-averaged, it is 2.5 for Egas < 10 Mqpc"^ and 
A'' ~ 1.4 for higher densities [20]. If disk thicknesses do not vary much among galaxies, 
then the SK law can be written in dcprojccted units with the same index A^. 

For a source with gas number density n, proton energy density Up, and volume V, 
the integrated hadronic 7-ray photon luminosity above some photon energy e is 

L>, = [ g>,nUpdV s'^ (2) 
Jv 

with the integral emissivity g>e in units of photon s~-'^[H-atom]~-'^[eV/cm^]~-'^ [9]. Using 
volume-averaged quantities and setting e = 100 MeV, from Eqs.(l),(2) we can write 

i>100McV -^'^gas^SN- (3) 

This is in agreement with the observational luminosity vs. (gas mass) x (SN rate) corre- 
lation [16], which evidently describes the 7-ray luminosity arising from tt*^ decay. 

Owing to the SK law, the above equation transforms intoL>e oc SFR^+i/^: if A^ = 2.5 
as appropriate for our sample galaxies with Egas < IOA/qpc"^ [20], then 

-^>iooMcV oc SFR^'"* . (4) 

Within limited statistics, this prediction agrees with observations [16]. 

So the observational non-linear L-^-SFR correlation [16] stems from the GeV lumi- 
nosity being (mostly) hadronic in origin, from CRs being linked with SF (through SN 
explosions), and from the Fermi/LAT-detected galaxies being located in the steep wing 
(i.e., low-Egas regime) of the SK law of SF. 
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5. Conclusion 

The link between SF and CR particles was suggested long ago [21]. It is based on the 
recognition that the CRe's, responsible for diffuse non-thermal synchrotron emission, are 
produced in the sites of SN explosions. The rough agreement between the relatively short 
lifetime (^S 3 x 10^ yr) of massive stars, and the similarly short synchrotron energy loss 
time of high-energy electrons, led to the expectation that the non-thermal radio emission 
of a galaxy is a measure of its SF activity on scales much shorter than the Hubble time. 

HE/VHE 7-ray detections (with, respectively, the orbiting Fermi telescope and ground- 
based lACTs) of some nearby star-forming (SB and normal) galaxies, that span a large 
range of SFR. have provided direct Up measurements. These are consistent with theo- 
retical predictions based on radio measurements, and with estimates based on SN rates 
and local CRp residence times. Should the match between measured and predicted Up 
be confirmed, some immediate implications would be: 

(i) star- forming galaxies can be powerful particle accelerators, able to achieve CRp energy 
densities orders of magnitude higher than the Galactic value; 

(ii) SNe, both in quietly star- forming galaxies and in very actively star-forming galaxies, 
probably have a common universal CR acceleration efficiency; 

(in) CR energy densities and equipartition magnetic fields derived from radio measure- 
ments can be used as proxies for the quantities characterizing the full particle energy 
distributions (derived from accurate spectral fits of the GeV-TeV emission): this could 
be particularly useful in the case of galaxies that are too far away for their (umbeamed) 
7-ray emission to be measured. 
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